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HEAT CONDUCTIVITY OF GRANULAR SYSTEMS 


by 

G. N. Dul'niev 
Z. V. Sigalova 


An analytical description of heat transfer in a model of a granular 
system is given. New expressions are derived for heat conductivity of granular 
systems with filling gas under normal and under various pressures. 


Most resembling the actual structure of a granular system would be 
a model with a chaotic distribution of particles having nearly equal volumes 
and rough surfaces. However, the difficulties connected with calculating the 
geometrical and physical parameters of such a system make it necessary to 
idealize the structure of the granular system and to assume for it the presence 
of a "distant" order in the distribution of the particles. Since the particles 
have a rough surface, 1 we will separate conditionally the spherical base in each 
particle and will regard the region occupied by the rough portions as an aureole 
of uniform thickness. 

Assuming that the most probable compact packing of the particles is 
tetrahedral, a sphere in such a system will rest on three spheres located below 
it. With a tetrahedral arrangement of the spheres, the porosity is equal to 
25. 95 percent and is independent of the size of the sphere [1] . In actual 
systems the larger porosity is due to the irregular shape of the particles, 
deviation from ideal packing, and other causes. 

Let us use a four-layer, most compact packing (the repetition of the 
first layer begins in the fourth layer) , and let us consider the transfer of heat 
in direction I (Figure la) . 

In any system with a long-range order we can treat separately the 
smallest volume (an elementary cell) whose innumerable repetition will give 
us the original system. 

In our system the elementary cell is represented by a rectangular 
parallelepiped ( Figure lb) whose base is a square with a side d = 2r and a 
height h - d*/2T 


1 


I 




FIGURE 1. MODEL OF A 
GRANULAR SYSTEM 

a) Most compact packing of spheres; 

b) elementary cell: 1) central solid 
particle; 2) eight equal parts of solid 
particle; 3) two halves; 4) four one- 
quarter parts of octahedrally-shaped 
large pores; 5) eight halves of 
octahedrally-shaped small pores. 



FIGURE 2. FOR CALCULATING 
THE HEAT CONDUCTIVITY 
BETWEEN SPHERES 


It is not difficult to prove that 
the effective heat conductivity of a 
system with a long-range order in 
distribution of the particles is 
identical to that of an elementary 
cell and, therefore, all further dis- 
cussions will deal only with 
elementary cells. 


We will show how the change in 
porosity in actual systems (p> 26 per- 
cent) depends on the size of the 
aureole a. ( Figure 2) . The volume 
occupied by solid particles in a 
granular system is equal to (100-p) 
and the volume of the gaseous aureoles 
is (p-26) . The ratio of the volume 
of the gaseous aureole to the volume 
of the solid particles expressed 
through geometric parameters is 
equal to : 


p-26 
100 - p 




Solving equation (1) for ( a/r ) , 
we obtain 



The following flows of heat 

( Figure lb) pass through an 

elementary cell by overcoming the 

resistance R.. Four identical flows 
i 

through a part of the solid particles. 
Four identical flows from the upper 
particles 2 to the central particle 1 
passing through: a) the gaseous 
aureole ( R t ) , b) the microclearance 
between the contacting rough portions 
(R 2 ) , and c) by direct contact with the 
particles R^. A flow to the central 


2 


particle 1 through one-half of the large pore 3 (R 3 ) . A flow through a part 
(Vg) of the lateral large pores 4 <R 4 ) . Since the resistance of the large pores 
is zero (an assumption) , therefore, all thermal resistances R. are parallel 
(the conductivities cr. = 1/R. are added together) . 


Assuming that the upper and central planes of the elementary cell are 
isothermic and taking into consideration that 4 <t 4 = a 3 , we obtain the following 
expression for the effective conductivity a of one-half of the elementary cell: 


a 


ecfr 

0.5h 


4A r A 

— 7 ^ 1 — = 4(0-! + (r 2 + 0. 5a 3 ) + 4<x . 

\T2 k 


(3) 


ecj) = effective 


It follows from equation ( 3) that 


% = 7a ,<r > + a ‘ * + V (4) 

e<p = effective 

Consequently, the problem is reduced to a determination of the 
structure of the heat conductivities. From now on, only dry granular systems 
will be considered. 

We will assume that the resistance of a solid particle is zero (the heat 
conductivity of the sphere's material is infinitely great) and that the entire 
temperature drop At takes place in the clearance. The thermal flux due to the 
molecular transfer through the clearance is equal to: 

dQi = ^ “ dS = X ~ xdxd<p . ( 5) 

r 2y r 2y 


r = gas 

Expressing y through x and integrating the equation ( 5) for the angle 
<p within 0 to 2t and for x from 0 to r, we can find the thermal flux through 
the aureole (the small pores in the cell) 

Qi = * r At7rr ^A In A “ 

r 
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On the other hand, we have Qj = Atov A comparison with (6) will give us 


a, = X f r Ain - 1 

r \ A - 1 


( 7 ) 


The heat conductivity a 2 of the gaseous clearance between the contacting 
rough portions can be expressed by the following relationship: 


<J 2 


A s u 
r H/ 


1. 5h 

ill 


H = initial 
Hi = sphere 
r = gas 


Let us first consider the contact of two ideal (free of microroughness) 
spheres with a radius r and let us find the drawing together of these spheres 
when acted upon by an outside load and forming thereby an area of a circle S. 
From the theory of elasticity it follows that the following relationships [2] are 
valid for contacting smooth spheres subjected to a load: 


9(1 - m 2 ) 2 P 2 


Vs 


b = 


H 


(1 - ff) 3P 2r 


-l 2 /s 


2rE 2 


S = 7T 


H 


8E 


( 8 ) 


H = initial 


The expressions ( 8) will have a different form when the centers of the 
spheres and the direction of the outside force P are not located on the same 
straight line. 

Assume that a certain load whose specific value is designated by A is 
applied normally to the upper plane of an elementary cell. The solid surfaces 
of the cell's upper plane have an area of 7rd 2 /4 and, therefore, the value of the 
load will be A7rd 2 / 4. This force is distributed among four particles. Taking 
into account the geometry of the elementary cell, the load per each contact 
will be 


P = A7Td 2 /l6^. 

H 

H = initial 


Let us assume that the elementary cell is not subjected to lateral 

deformation, i. e. , p = 0, and that the force P produces only a longitudinal 

H 

deformation; in this case we obtain from equation ( 8) 
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b = 0. 56d 


(9) 



S = 


0. 14 7rd 2 



Let us introduce the concept of the modulus of elasticity E 0 of the 
granular system consisting of compactly packed smooth spheres. The relative 
deformation D of the elementary cell acted upon by the specific load A is equal 
to 


D 



2b 
h ' 


From the equations (9) and (10) we obtain the relationship 


( 10 ) 


V 2 / 

E n = 1. 3 A /3 E . 


Experiments had shown that 


E 0exp = * - 0.3 f 0.6, 


(ID 


i.e. , $ + 1. 3, as it follows from the theory. The possible reasons for this 
discrepancy are as follows: The particles are not arranged as required by an 
ideal theoretical pattern; the particles have a wavy and rough surface. With 
this in mind, we will use the expression (11) for approximated calculations. 


Taking the relationship (11) into account, we obtain 

4/ 


S H - 


0. 14 7r d^ 


A* 


*) • 


H = initial 


Since the area in actual contact represents an insignificant part of the nominal 
area (10 -2 to 10 -5 ) , we will assume that the area of the gaseous layer is equal 
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If the material of the granular system received a preliminary treatment, 
the height of the sphere h^will be determined by the class of cleanness of its 

surface. If the surface of the material was not first treated, the height of the 
sphere h^ is determined experimentally. Tests made with gravel and lead-shot 

make it possible to expect a stable h^/d ratio, namely, 


h m /d = y • 10 — 3 ; y = 1 -r 3. 

i ii = sphere 


Substituting the values of h and S in the value of <r 2 we obtain an 

ILL H 

expression for the conductivity of the gaseous clearance 


<?2 


, ^d . io 2 

r 2 4 

r T/T 3 



r = gas 


( 12 ) 


The heat transfer through the direct contact X ( through the hard body 

xv 

of the micro roughness) is still determined only experimentally. An analysis 
of the experimental investigations shows that X = 1 x 10 -3 to 5 x 10" 3 for 

K 

nonmetallic particles and 0. 3 to 0. 5 for metals. 


Into the large pore can be inscribed a sphere with a radius of 0. 41 r [1] . 
Taking into account the gaseous aureoles, the radius / r \ of the inscribed 
sphere will be equal to ' 

r = 0. 41 (r + a) + a. 

o 


Conditionally, the large pore will be regarded as a cylinder with a 
radius of r and a height of 2r . The molecular conductivity of such a cylinder 
is equal to £T 

<r 3 = X 2 7ir(l. 41 A - 1) . 
r 

With this expression the conductivity a 3 is calculated for a porosity when the 
radius of the pore becomes equal to the radius of the particle, which 
corresponds to a/r = 0. 41 or to p = 75 percent. 
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For a porosity larger than 75 percent there appears a new channel of 
the transfer: a column of gas with a height of rAN/iTand an area of 2r 2 (2A 2 - it) . 
Its molecular conductivity <t 5 is equal to 


^5 


2r(2A 2 - 7 r) 

r 4Ya 


(13) 


Substituting in equation (4) the values of a. from equations (7) , (12) , 
(13) , and it is possible to obtain an expression for the effective heat con- 
ductivity of the granular system. This expression, however, will still fail to 
reflect the heat transfer through the system because, as above, a number of 
essential phenomena have not been taken into account. We will proceed with 
the analysis of these phenomena. 

a) The calculation of the conductivity of the gas cr 2 between the 
microrough portions did not take into account the dependence of cr 2 on the 
porosity, although it is obvious that an increase in porosity must reduce the 
value of u 2 . Let us use cr 2 for the true value of the sought conductivity, in which 
case 


cr 2 = cr 2 f (p) . (14) 

It is obvious that f(26 percent) = 1 and that f(100 percent) = 0. 

Let us present f(p) in form of a difference between two functions 
26/p - (p(p) and let us assume that cp(p) varies linearly with p when passing 
through two points (p( 26) = Oandcp(lOO) = 0.26, in which case 


f(p) 


— - 0. 0035p + 0. 09. 
P 


b) Let us estimate the amount of heat carried through the gas by 
radiation. As shown by A. F. Chudnovskiy, the radiation component X of an 

th P 

effective coefficient of heat conductivity in the i clearance can be calculated 
by using the expression [3] : 


X . = 2e 2 CT 3 6 . = gS ., 
pi p mi mi 



(15) 
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th 

The additional conductivity a . due to the radiation in the i gaseous 
element of the cell is equal to P1 

a . = X .S ./ 6 .. 
pi pi mi mi 

In calculating the effective heat conductivity of the system by using a 
previous formula [4], it is necessary to add cr . to <j. for i = 1, 2, 3, and 5. 

pi i 

c) The molecular heat transfer through a gaseous layer depends on 
many factors. As has been shown [ 4] , the heat conductivity of a gas ^ X^ whose 

pressure varies widely can be calculated with the aid of the formula 


X = X n 


1 + 5 i rT< Pr >'' 


M 

r = gas 


Kn 


-l 


(16) 


H, = 


If A t is the length of the free path of a gas molecule under a pressure of 
133,322 H/m 2 , then 


A = 


_ AlH! 


H 


(17) 


Taking into account the relationships (16) and (17) , the expression for 

.th 


the heat conductivity X in the i layer can be written as: 


X n 


X ri 1 + B/H5 . 

mi 


, B = —^-(Pr)- 1 (^-^) AiHt. 


k + 1 


a 


It can be shown that in an elementary cell the expression for a . and 

S . will be as follows: mi 

mi 


ml 


5 _ = h , S 0 
m2 ill m2 


£ 

4 


’ S ml 

7rr 2 ; 

0. 14 7rd 2 

(e\\ 


V E / • 
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6 m3 


S 

2r(l. 41A - 1), g 

6 c = 2 \HfrA, 
m5 


m3 

m3 

S m5 


7i r 2 ( 1 . 41A - l) 2 , 
7i r 2 , p & 75%; 


= 2r 2 (2A 2 - 77). 


P =£ 


75%; 


ill = sphere 

The final expression for an effective heat conductivity of a granular 
system assumes the following form: 

at 26% < p < 75% 

X (x + Y + 2.23 1 14 ^ h ~ 1 ) 

e0 A y 1 + B/HS^g ) 


4 45 n 

+ gr [1 + 0. 5(1. 41A - l) 2 ] + \ K ; 


at p > 75% 


(18) 


= ^- j^X + Y+2.23 


e<£ A 


(1. 41 A - 1 ) (1 + B/H6 m3 ) 


where 


A 2 - 1.58 


] 


4. 45 (. AM 

(1 + B/H • 2. 82rA) A J + A gT \ + n J + 

ec p = effective 


X = 


1 + B/Ho 
Y = f(p) 


4. 45 A 

x Ain 


ml 

1 


A - 1 


- i ; 


9.9 • 10 2 /a\^ 


1 + B/H6 


m2 


T/3 


% 


The formulae expressions (18) were used to calculate the coefficients 
of heat conductivity of various granular systems for a wide range of porosities 
and pressures of the filling gas. A satisfactory agreement was obtained 
between the calculated and experimentally checked values of the coefficient of 
heat conductivity of granular system ( Figure 3) . 
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DESIGNATIONS 



FIGURE 3. COMPARISON 
OF EXPERIMENTAL 

(1.3) AND CALCULATED 

(2.4) DATA 

1 ) steel balls with a diameter of 
1. 26 mm in hydrogen at p = 35 per- 
cent [5]; 3) silica gel, d = 5-6 microns 
in air, p = 95 percent [5] ; 2 and 
4) steel balls and silica gel. 


d — diameter of a particle; 
h — the height of elementary cell; 
2 a; — the aureole thickness; 
p — the porosity of the system; 

2 y — the distance between 
particles; X e( p~ the effective heat 

conductivity of an elementary cell; 
X — the molecular heat con- 


r 

ductivity of the gas in the clearance; 
Xq — the heat conductivity of the 
gas under normal pressure; 

R. and a. — respectively, the heat 

resistance and conductivity of the 


i part of the elementary cell; 

cr t and 05 — the conductivity of the 

gaseous aureole and of the large 

pores; a . — the average size of 
mi 

the clearance, in meters; 

S . — the average area of the 
mi 

i^h element; S — the nominal 
H 


area of the contact (the area of the gas layer between the contacts) ; h^ — the 


average height of the clearance (between the microrough portions) ; b — the 

convergence of the spheres; P — the contacting force; p — Poisson’s ratio of 

H 


the material of the particles; E — Young's modulus of the material of the 
particles; e — degree of blackness of the surface of the pores; C — a Stephan- 


Boltzmann constant; T — the average absolute temperature of the material; 
k = c / c^ — the ratio of heat capacities of the gas at constant pressure and 

volume; a — coefficient of accommodation; A — the length of the free path of the 
molecules; H — the pressure of the filling gas. 
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